Multiple scattering of light in cold atomic clouds with a magnetic field 
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Starting from a microscopic theory for atomic scatterers, we describe the scattering of light by a 
single atom and study the coherent propagation of light in a cold atomic cloud in the presence of 
a magnetic field B in the mesoscopic regime. Non-pertubative expressions in B are given for the 
magneto-optical effects and optical anisotropy. We then consider the multiple scattering regime and 
address the fate of the coherent backscattering (CBS) effect. We show that, for atoms with nonzero 
spin in their ground state, the CBS interference contrast can be increased compared to its value 
when B — 0, a, result at variance with classical samples. We validate our theoretical results by a 
quantitative comparison with experimental data. 
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I. PHYSICAL CONTEXT 

Multiple scattering of waves is an important topic in- 
volved in many branches of physics, from atomic physics 
to astronomy via condensed matter physics [2 [5] . One 
of its most fascinating aspects occurs when interfer- 
ence effects persist in the presence of disorder. These 
interference effects manifest themselves through devia- 
tions from the usual diffusive behaviour obtained at long 
times (weak localization) . Under suitable conditions they 
can even completely inhibit transport (strong localiza- 
tion) reaching an insulating regime. During the last two 
decades, both the weak and strong localization regimes 
have been extensively studied with systems ranging from 
electronic waves to ultracold atoms [3til2). 

Coherent multiple scattering of light in cold atomic 
clouds has been experimentally and theoretically studied 
since ten years in the context of the coherent backscat- 
tering (CBS) effect p3]. Compared to classical scattering 
media, such as semi-conducting powders, cold atoms con- 
fined in a magneto-optical trap (MOT) constitute quite 
a unique sample of identical strongly-resonant point-like 
scatterers [14]. A key feature of the atomic scatterers 
is the presence of quantum internal degrees of freedom 
when their ground state is degenerate i.e. possesses a 
non-zero angular momentum |15) . These internal degrees 
of freedom are coupled to the polarization of light waves 
during single scattering events. As they are traced out 
when one observes the interference pattern of light, light 
experiences decoherence during its diffusive propagation 
in the atomic cloud. This leads to a reduction of the in- 
terference between scattering amplitudes as exemplified 
by the CBS contrast [TB] . The saturation of the atomic 
transition |17. and the finite temperature of the atomic 
gas [Hj are two additional mechanisms reducing the co- 
herence of light. To restore the contrast of the interfer- 



ences, it has been suggested to use a polarized atomic 
gas [19]. Unfortunately, this method is efficient only in 
optically thin clouds. Another possibility is to break the 
degeneracy of the atomic ground state by applying an 
external magnetic field B, which is one of the topics we 
address here. 

In classical scattering media, the presence of an exter- 
nal magnetic field is known to modify the interference 
pattern and to reduce the CBS interference contrast [20] . 
The mechanism at work here is the modification of the 
polarization of light during its coherent propagation due 
to the Faraday effect [U [22] . As this modification is not 
the same along a scattering path and its reversed part- 
ner, a dephasing process takes place between the scatter- 
ing amplitudes associated to these two paths, scrambling 
in turn their interference. Although the magneto-optical 
effects are very strong in cold atomic gases [53], they 
actually have a negligible impact on the shape of the in- 
terference pattern itself [24j. More surprisingly, a full 
restoration of the interference contrast is even possible 
under well chosen conditions |25j . In this article, we give 
a detailed explanation of this result. We present a gen- 
eral, non-perturbative study of the interaction between 
a cold atomic cloud and quasi-resonant light in the pres- 
ence of a magnetic field. 

The paper is organized as follows. In section [llj we 
derive the scattering operator of a single atom exposed 
to monochromatic quasi-resonant light in the presence of 
a magnetic field B and we analyze the differential and 
total cross-section. In section IIIIl we address the coher- 



ent propagation of light and derive the refraction index 
tensor of the effective medium. In particular, birefrin- 



gence and magneto-optical effects are studied in HID 



and IIIF The CBS interference effect under a magnetic 
field is discussed in section IIVI and the restoration of the 
CBS contrast is studied in section fVl 
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II. SCATTERING OF LIGHT BY A SINGLE 
ATOM IN THE PRESENCE OF A UNIFORM 
MAGNETIC FIELD 

A. Physical setting and basic assumptions 

We consider the situation where a gas of identical 
atoms, with spatial density n, is both exposed to a uni- 
form external magnetic field B and a monochromatic 
light plane wave with wavelength A, wavevector k (k = 
27r/A), angular frequency w = k (we use units such that 
h = c — 1) and transverse polarization e. 

We first assume that the average distance between 
atoms is much larger than the wavelength A of the light 
wave, nX^ <C 1, meaning that multiple scattering pro- 
cesses take place in the dilute regime. In this case a 
semi-classical description is appropriate and propagation 
of light in the atomic medium is well described by par- 
tial waves propagating along classical multiple scattering 
paths. Each path consists in a succession of independent 
scattering events by a single atom, separated by propa- 
gation in an effective medium with refractive index Nr. 

To further properly describe scattering of light by one 
single atom, we assume that uj is very close to the tran- 
sition frequency ujq between the atom internal ground- 
state (with total angular momentum F) and an inter- 
nal excited state (with angular momentum F^), hereby 
considering the case of quasi-resonant scattering. Intro- 
ducing the light detuning S — co — loq, this condition 
reads \S\ <^ u/q. We further assume that this transition is 
closed (the excited state with life time l/F can only de- 
cay by spontaneous emission into the groundstate) and 
well isolated from any other allowed transition connect- 
ing the groundstate to any other hyperfine levels. This 
is experimentally achieved for instance in the case of the 
F = 3^ Fe^A transition of the D2 line of ^^Rb [13 . 

The total atom-light Hamiltonian H = Hq + U is the 
sum of the free atom-light Hamiltonian Hq and of the 
interaction term U describing the coupling of the light 
modes to the atomic degrees of freedom. When no mag- 
netic field is applied, the interaction between light and 
one single atom is well documented [2^1 HZ] and we thus 
simply need to incorporate the effect of the magnetic 
field. We will assume here that weak fields are applied 
to the gas so that the linear Zeeman effect is the relevant 
physical description. As the magnetic field only affects 
the energy levels of the atom, the free atom-light Hamil- 
tonian Hq just reads: 

Hq = ujQpe + ^.gB ■ F + figeB • Fe -f ^ i^kai^ake- (1) 

fc.e 

Here /j, « 1.4 MHz/Gauss is the Bohr magneton and g 
and ge are the Lande factors of the ground and excited 
states respectively. In the following we will choose B to 
be along axis Oz which we choose as the quantification 
axis for the internal Zeeman states of the atom. Then 
Pe = J^m =-F l-^eTOe) (-FeWel IS the projector outo the 



excited state where \Feme) denotes an atomic state of 
angular momentum Fg and magnetic quantum number 
nie referred to axis Oz. 

To be consistent with our approximations, the Zeeman 
shifts (of order fiB) must be much smaller than the en- 
ergy difference between any other atomic hyperfine levels. 
Still these shifts can be sufficiently large to fully split the 
Zeeman structure and strongly modify the light scatter- 
ing properties of the atom. Indeed, as exemplified by 
the case of Rubidium atoms, Zeeman shifts comparable 
to the excited-state frequency width {T/2tt — 5.9 MHz), 
are achieved with moderate field strengths of 5 Gauss. 

As usual flfcg and a]^^ in (llj) are the annihilation 
and creation operators of the electric field mode with 
wavevector fc, polarization e _L fc and frequency ujk. We 
further assume that the light wavelength is much larger 
than the size of the relevant atom electronic wavefunc- 
tions. The interaction between light and an atom located 
at position r is then accounted for in the dipolar approx- 
imation and reads U{r) = — dg • E±{r). The atomic 
electric dipole operator writes = d + d) where 

= P,d,Pg (2) 

is the operator describing transitions from the ground 
state to the excited state. Here Pg = X]m=-_F 1^™) {Ffn\ 
is the projector onto the groundstate sector with \Fm) 
an atomic state of angular momentum F and magnetic 
quantum number m referred to axis Oz. The transverse 
radiation field operator reads E^_{r) — D{r) + D^{r) 
where: 

£)(r)=j^£^,efeafc,e*'- (3) 

k,e 

describes photon annihilation in all possible field modes. 
The field strength is S^^^ = yJuj^J^eoV , where V is the 
quantization box volume (it disappears at the end of all 
physical calculations). 

The dipolar interaction operator \J describes the pos- 
sibility for the atom to absorb or emit a photon changing 
at the same time its internal state \Fm) — \F(m^ or 
\Ff,me) — \Frfi). As the incident light is nearly resonant 
with the atomic transition, we only need to consider res- 
onant contributions where the atom in its ground state 
can only absorb a photon, and the atom in its excited 
state can only emit a photon: 

U(r)KD{r)-S ^D\r)-d. (4) 

This is known as the rotating wave approximation. 

It has to be noticed that the photons considered in this 
article are associated to the transverse electric field, i.e. 
to the electric displacement vector D and not to the elec- 
tric field E (hence our notation) This point is im- 
portant for two reasons: first, due to the magnetic field, 
light propagates in an anisotropic medium where D and 
E are no longer colinear and, second, former theoretical 
studies treated the propagation of E |22J. 
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B. Scattering amplitude 

Starting from the incident state \i) — \kijje; Fm) , de- 
scribing a photon \kuje) impinging an atom located at r 
with initial ground state \Fm), we consider the scatter- 
ing process leading to the final state |/) = \k'uj'e' ; Fm') , 
describing a scattered photon \k'uj'e') leaving the atom 
in the final ground state \Fm'). The probability ampli- 
tude for such a transition is the matrix element (/| S \i) 
of the scattering operator 5 = 1— 2z7rT acting on the 
atom- photon Hilbert space H = Hat ®T~Ll- The transi- 
tion operator for an atomic point-dipole scatterer writes 



matrix element then reads: 



T : 



2V 



\r){r\ (g) T 



(5) 



where T couples the photon polarization to the atomic 
internal degrees of freedom. The (on-shell) matrix ele- 
ments of T are: 



i{k—k')-r 



6{E' - E) Tf,{E) (6) 



where the delta distribution ensures energy conserva- 
tion. For the initial and final states under consideration, 
E' = Lo' + iigBm' and E = uj + j.igBm. It is important to 
note at this point that, contrary to the case where there 
is no magnetic field, scattering has now become inelastic, 
a feature that will complicate greatly the analysis of the 
multiple scattering situation. Indeed as soon as m' ^ to, 
the atom experiences a net Zeeman energy change. In 
turn, the angular frequency of the scattered photon is 
accordingly modified to w' = u) + fig B {m — m') . It should 
be noted that we simplified here the problem further by 
neglecting momentum transfer to the atom during scat- 
tering. This approximation, supposedly valid when the 
recoil energy is negligible compared to the energy width 
of the excited state, becomes nevertheless questionable 
when considering the multiple scattering regime as stud- 
ied (without magnetic field) in [H]. 

The matrix element 7/i {E) is calculated from 7/i {E) = 
(/I UGe{E)U\i) where Ge{E) = P^{E - H)-^Pe is the 
Green's function of the system when the atom is in its 
excited state. This operator can be computed by using 
resolvent techniques [27] and we find: 



Ge{E) = 



E 



|FeTOe) (FeTOe 



E — geTTlellB — (jJq + iT /2 



(7) 



where F is the angular width of the atomic excited state 
due to coupling to vacuum fluctuations. Technically 
speaking, there is also a modification of the atomic fre- 
quency (Lamb-shift) that we get rid off by a proper re- 
definition of ujq . It is noteworthy that neither the Lamb- 
shift nor the linewidth depend on B. This is so because 
the Zeeman operators do not couple to the dipole inter- 
action U. 

Introducing the reduced atomic dipole operator dg — 
de/de with dgWQ = STreoc^^F |27| . the internal transition 



r (F)-^"" r/2 



m' 7n 7 



(8) 



where the dyadic transition operator in polarization 
space tjn'm is: 



tr 



^ {Fm'\ de |FeTOe) (-FeWel dg \Fm) 



The dimensionless parameter 

, 4'B 



1 - i(t){gm - geirie) 



1 - 2iS/T' 



where 



2/iS 

1^' 



(10) 



(11) 



quantifies the impact of B on the atomic scattering prop- 
erties and is generally complex-valued except at reso- 
nance (<5 = 0). One can note in passing that, when B — 0, 
tm'm — {Fm'\dS \Fm) which describes the usual ab- 
sorption and emission cycle from the degenerate ground- 
state [551. 



C. Scattering cross section 

Probability conservation assures that, under the scat- 
tering process, any incident photon is either transmitted 
in the same mode or scattered in another mode. This 
is the essence of the optical theorem, which relates the 
total scattering cross section of the atom to the forward 
scattering amplitude: 



cr = -2VIm {kuje\ (r),^„j \kue) 



(12) 



where — Tr(» pat) indicates an average over the 

initial atomic internal degrees of freedom. As such 
2V (r)j„t = <8) (T)j„f only acts on the photon de- 

grees on freedom. At this point, we make the impor- 
tant simplifying assumption that the initial atomic den- 
sity operator pat describes a complete incoherent mix- 
ture of the groundstate Zeeman sub-states. This is gen- 
erally the situation for cold atoms prepared in an opti- 
cally thick magneto-optical trap (MOT) where all Zee- 
man states are uniformly populated and no Zeeman co- 
herence is achieved. The initial atomic density operator 
then reads pat — Pg/C^F + 1)- Using the fact that de, e 
and e' are irreducible tensors of rank 1, one can express 
e' • (i) • e in terms of the irreducible components of the 
tensor e'i£j in the cartesian basis set {x, y, z). Introduc- 
ing the dyadic projector = 1 — aa/a" onto the plane 
perpendicular to a, we find: 



M, 



2K 



Gtt F/2 
S + ir/2 

F 

m=-F 



FF^ 



(13) 



(14) 
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where MpF, = (2Fe + 1)/(3(2F + 1)) is a factor tak- 
ing care of the degeneracies of the ground and excited 
states. Both projectors Aj- and A^/ ensure that polar- 
ization vectors always remain transverse to the direction 
of propagation. 

Using the totally antisymmetric tensor of rank 3 Sijk 
{sxyz = 1), the dyadic tensor T"*^^^"' reads: 



observed when k B 



rjSijkBk + S,BiBj. 



(15) 



With the z-axis chosen along B, the matrix representing 




(16) 



The Cj V ^nd ^ coefficients depend on (f> and on the 
Clebsch-Gordan coefficients of the atomic transition 
HO]: 



lit, 

\ m=-F 



+ 



{FF, 



mm 



'.)m- ge) 
-111-1)2 



n=-F 



1 - i(t){{g - ge)m + ge) 
(FFe- mm -hl|ll)2 



i<f>{i9 - 9e)m - ge) 



e = -c- 



{FFe - mm - l\l - 1)^ 
1 - 'i<l>{{9 - 9e)m + 9e) 
^ {FFe-mm\W)'^ 

m=-F 



i(l){g - ge)m 



(17) 



As is expected from Onsager's reciprocity relations [5T] . 
one can check that ( and ^ are even function of (j), while 
r] is an odd function of (j). It can be checked that all 
these three coefficients are real at resonance {6 = 0) since 
(j) = 4>B then real. In the case ofaF = 0— >Fe = l 
transition (17 = 1), these coefficients read 



c 



»7 = 



1 



(ffe0)2 

9e(t) 
+ (5e0)2 



(18) 



1 + {9e4>)^ 



The dyadic tensor T'^^^'^ embodies the effect of the 
magnetic field on the photon polarization degrees of 
freedom and gives rise to the usual magneto-optical 
effects. The C term is responsible for normal extinc- 
tion (Lambert-Beer law). The 77 term describes the 
magnetically-induced rotation of the atomic dipole 
moment (Hanle effect) [211 [32] and induces Faraday 
rotation and dichroism effects observed when fe || S [33]. 
The ^ term is responsible for the Cotton-Mouton effect 



From expression (13) and the optical theorem (12) we 



deduce the total scattering cross section of a photon by 
an atom initially prepared in an incoherent mixture of 
Zeeman internal ground states: 



a((/)) = a Re (^(1 + 2^^) (e • T^^^^^ ■ e)^ 



(19) 



(J Re ( (l-f 2^-)(C + rj(e x e) ■ B + £,\e ■ Bp 



where 



In the case ofaF = 0— >Fe = l transition {g — 1), and 
for resonant light {5 = 0), the scattering cross section 
boils down to the simpler form: 



(j{cj)) = (To 



l + {ge<t>B)^\e-B\ 

{9e(l)B)'^ 



1 



(21) 



As one can see, in the presence of a magnetic field, the 
total scattering cross section depends explicitly on the 
incident polarization e. More precisely, it depends on 
the relative direction of e with respect to B. We recover 
here the well-known fact that an external magnetic field 
induces optical anisotropy in otherwise isotropic media. 
In the absence of a magnetic field, T'^^^''' reduces to the 
identity matrix 1 and we get (t(0 = 0) = tr, giving back 
the result found in [2SI- 



D. Impact of optical pumping 

Our previous calculation in fact just considered the 
scattering of one quasi-resonant photon by one atom as- 
sumed to be initially prepared in an incoherent mixture 
of Zeeman internal ground states. In a real experiment 
however, many quasi-resonant photons are shone. For 
a given incident polarization, the repeated scattering of 
photons by the atom induces changes in the populations 
of the groundstate Zeeman sublevels and creates coher- 
ences between them. This effect is known as optical 
pumping and is enhanced in the presence of a magnetic 
field. Our previous calculations thus applies to the case 
where optical pumping can be neglected. This is the sit- 
uation considered in section |Vl One can note however 
that our results can be easily extended to include optical 
pumping. Indeed, the most general tensor of rank 2 that 
can be written with the components of B is still given by 
expression ( 15 ) but with arbitrary coefficients now. Thus 
equations (13), ( [T6| an d (19) remain valid. Only the de- 
ms (|17[) of C, »7 and f will be modified. In 



tailed expressions \Vl\ of 77 and ^ will be modified, 
the following, all our analytical results are expressed in 
terms of these three coefficients. They thus remain valid 
in the presence of optical pumping if the coefficients are 
given their appropriate expressions and values. 
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E. Differential scattering cross section 

We now turn to the impact of B on the radiation pat- 
tern of the atom. The intensity of hght scattered in the 
direction k' with polarization e' while the atom changes 
its internal state from \Fm) to \Fm') is proportional to 
the differential scattering cross section 



(27r)2 
9a 2F+1 
8^ 2Fe + 1 



k'e',Fm'\T\ke,Fm) 



e' ■ t„ 



(22) 



Starting with an atom prepared in state \Fm) the differ- 
ential cross-section for a photon to be scattered in mode 



dfl 



ike^k'e')= 



dn 



-(fee 



(23) 



The total photon scattering cross-section is then obtained 
by averaging over the initial internal atomic density ma- 
trix, which we assumed to describe a fully incoherent 
mixture of groundstate Zeeman sublevels. We thus ar- 
rive at: 



da 

dn 



(fee ^ fe'e') 



1 



2F 



9a 



F 

m.m' — — F 



dan 



dVL 



(ke k'e' 



E 



^71 



2Fe + l 



(24) 



When the ground state is non-degenerate {F = 0, i^e = 
1, g — 1), the tensors ^oo and T^^^-* coincide and the 
differential cross section takes the simple form: 



da 

dn 



(fee ^ k'e') 



^1^7.7^(01) 

Stt 



(25) 



medium by a homogeneous effective medium having a 
complex refractive index tensor Nr. When B = 0, Nr is 
a scalar and bears no action on the incident polarization. 
Its imaginary part gives rise to an exponential attenua- 
tion (Lambert-Beer law), its characteristic length scale 
being known as the extinction length. Since true ab- 
sorption (i. e. conversion of electromagnetic energy into 
another form of energy) is absent in our case, depletion 
of the incident mode can only occur through scattering. 
The extinction length thus identifies with the scattering 
mean free path £ and one has Im(iVr) = l/{2k£). 

The presence of a magnetic field B modifies the inter- 
action between light and atoms and the coherent prop- 
agation of light in the atomic cloud will be accord- 
ingly altered. As B induces a preferential orientation 
of space, it will set an optical anisotropy in the atomic 
gas. In this case, the refractive index tensor is no longer 
a scalar and will act on the polarization space: the co- 
herent propagation of light will exhibit birefringence and 
magneto-optical effects. These magneto-optical effects 
are well documented in the literature [2TH23] : for a given 
wavevector fc, there are two eigen-polarization modes 
propagating in (possibly) different directions, with dif- 
ferent velocities and attenuations. 

In this paragraph, we extend the techniques developed 
in and [51] to the case of atoms with a degenerate 
ground state. We will always assume that the atomic gas 
is dilute. In this case, nX'^ <C I (n being the number den- 
sity), and £ — l/{na). Furthermore, since for resonant 
scatterers a is at most of the order of A'^, £ will be al- 
ways much larger than the average interatomic distance 
n~^/^. As a whole, for point-dipole resonant scatterers, 
the dilute medium condition nX^ ^ I implies the con- 
dition k£ 1. The properties of the effective medium 
will then be directly related to the individual scattering 
properties of the atoms under the magnetic field. 



with 



+ ry (e' X e) 5 + ^ (e'- B)[e- B), (26) A. Average Green's function for light propagation 



the coefficients (, t] and ^ being given by ( [l8| . 

Since the resonant denominator in i|9l explicitly 
depends on the magnetic quantum numbers m and mg, 
standard irreducible tensor methods [3S] are of little 
practical use to boil down the total differential cross- 
section ( 24 ) to a much simpler form. However it can be 



easily and efficiently computed for any transition line 
by using the symbolic calculus software Maple"^^. As a 
humorous note, the reader is invited to appreciate the 
power of the optical theorem by deriving the total scat- 
tering cross section ( 19 ) by direct computation from ( 24 1 . 



III. COHERENT PROPAGATION OF LIGHT 

The propagation of an incident light field mode \ke) in 
a scattering medium, also known as the coherent propa- 
gation, is efficiently described by replacing the scattering 



The first step to find the refractive index tensor Nr and 
describe the coherent propagation of light is to determine 
the Green's function for light once all atomic degrees of 
freedom have been averaged out. 

Let N be the total number of atoms in the gas and 
their respective positions be labeled by {i — 1,- ■ -N). 
The Hamiltonian of the system {atoms -I- light} is 

N N 
i=l 1=1 

= Ho+U, (27) 

where i?o a-nd U are given by expressions (fTl) and Q . The 
Green's functions of the whole system G(z) = [z — 'H)~^ 
and of the uncoupled system Go{z) = [z — Ho)^^ satisfy 
the recursive equation 

G{z) = Gq{z) + Go{z)UG{z) 

= Go{z) + Go{z)UGo{z) + . . . (28) 
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Since we are interested in the situation where all atoms 
start and end in their groundstate, we merely look for 
the Green's function projected onto the atomic ground- 
state manifold PgG{z)Pg. In this case, only expansion 
terms containing lA an even number of times can con- 
tribute. The average over the atomic degrees of free- 
dom will generate a Dyson equation for the average 



modes are: 



Green's function Q{z) 



{PgG{z)Pg). Under the di- 



lute medium assumption, and since all atoms are iden- 
tical and uniformly distributed in space (at least on the 
scale of the scattering mean free path), the correspond- 
ing self-energy is given by = N{Ti{z))int where 
Ti{z) = U{ri) + U{ri)Go{z)U{ri) -t- ... is the transition 
operator of the i-th atom. The average over the internal 
degrees of freedom is given by (13 1. 



The matrix elements of the average photon Green's 
function are then given by: 

{k'e'\giu;)\ke) = 6kk' ^ ■ Q{k,Lo) ■ e (29) 

where — l/(n(To)- The Kronecker symbol 5kk' fea- 
tures the restoration of translation invariance under the 
spatial average. The self-energy tensor S(/c, k) contains 
all the information on the effective medium. It has an 
explicit dependence on the incident direction because, in 
the presence of the magnetic field, the scattering medium 
develops an optical anisotropy (see section II C ) . 



B. Optical anisotropy 

The equation ( [30j ) receives a simple interpretation in 
a basis where I](fc, fc) is diagonal. It corresponds to po- 
larization modes which propagate in direction k without 
deformation. The poles of Qik.uj) then give the corre- 
sponding dispersion relation for the eigenmode. 

Trivially, k is an eigenvector of E(A:, fc) with eigenvalue 
0. This is a consequence of the transversality of light, 
and this eigenvector is not physically relevant. To find 
the other complex eigenmodes and complex eigenvalues, 
E(/c,fe)V± = A-t'Vj-, we parametrize k by its spherical 
angles and in a coordinate frame with z-axis parallel 
to B. We find: 



A±(fc) = C + e 



sm 



± \l~-n'^ cos2 6* -h ^2 



sm 



(32) 



Their explicit dependance on the angle 9 between B and 
fc expresses the optical anisotropy of the atomic cloud 
induced by the magnetic field. The corresponding eigen- 



V± (k) cx cos^ 9 cos ip + £, sin f 



=F sin (p\ —rf cos^ -I- 



/ • .sin^e 

- ^77 cos d sm — ^ — ^ — cos if 



, sin-* 9 



± cosi^ay —if cos^ I 
-?7cos0sin 9z 



,sin^0 



y 

(33) 



To allow a relative ease of reading, these vectors have not 
been normalized. 

One difficulty arises when {—rf cos^ 9 + £f{sm^ 9) /4) = 
for non- vanishing 77 and ^. This situation can only hap- 
pen if rj and ^ are real, thus for </> real, i. e. at resonance 
{5 — 0). The solutions are of the form {±9^,7^ ± 9o), 
which means, because of the invariance around the z- 
axis, that the photon has to propagate along a cone of 
apex 6*0. When this is the case, A+ = A_ but = VL, 
and the self-energy is non-diagonalizable. However this 
situation is largely unphysical in the sense that it arises 
from the first-order approximation in the atomic density 
used to compute Ti{k,k). At next order in density, this 
difficulty disappears. However, this could lead to inter- 
sting effects for the propagation of light near the apex 
angle. We chose to neglect such effects in the following. 

Noticeably, V+ and V- are not orthogonal vectors in 
general. Indeed, scattering depletes the coherent mode 
and its energy decays. This is reflected by the fact that 
the self-energy is not a hermitian operator. Thus its 
eigenvectors have no reason to be orthogonal. When 
-8 = 0, T^^^^^ reduces to the identity and E(fc, fc) is then 
proportional to the projector A^. In this case alone, all 
polarization states orthogonal to fc are eigenmodes and 
it is then possible to choose an orthogonal basis. One 
can however note that V+ and V are nearly orthogo- 
nal when \ri\ <C |^| or when \r]\ 3> |^|. This happens in 
the limit of small or strong magnetic fields {fiB <C S, F, 
fiB ^ S,T), or at a very large detuning {6 ^ ^B, T). 



C. Refraction index 

In the polarization eigenbasis, the poles of the Green's 
function (30) give the dispersion relation for V±. We 
find: 



uj± (fc) ~ k + 



1 



r/2 



2£o 6 + iT/ 2 



A±(fc) 



(34) 



The refraction index tensor Nr is diagonal in the polar- 
ization eigenbasis and the polarization vectors V± prop- 
agate each with different complex refractive indexes: 



uj±{k) 



1 



1 



r/2 



2k£oS + ir/2 



A±(fe) (35) 
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since Mq ^1. As a consequence, the two eigen- 
polarizations propagate with different phase velocities 
and experience different attenuations (dichroism) . fn 
turn, the effective medium acts as an absorbing polar- 
ization filter for the incoming light. 

The index mismatch between the two eigen- 
polarizations is 

^ (A+-A_)(36) 



2k£o 1 - 2id/T 




—rf cos^ d + 



sm 



(37) 



2Mq 1 - 

and vanishes for 9 = 6^, i. e. when the two eigen- 
polarizations collapse onto each other, rendering the re- 
fractive index tensor no longer diagonalizable. 

D. Group velocity and birefringence 



The Iui{r) integral has an ultra-violet divergence (fc — )■ 
oo) which needs to be regularized. It expresses that the 
interaction between two atoms separated by less than one 
optical wavelength cannot be reduced to the exchange of 
one resonant photon. The divergent part of the integral 
comes from the contact term of the radiated field |26] . As 
the average interatomic distance is much larger than the 
optical wavelength, we are interested only in the far-field 
component. The latter can be obtained from the regular 
part of ( 39 1 and calculated with the help of the residue 
theorem. The final result reads: 



Gu^[r) ^ -—Are 
Zirr 



ikrN, 



featuring the refractive index tensor 



Nr{k,f) 



A, 



Ar 



(41) 



(42) 



Let us consider a polarized monochromatic wave 
packet propagating in the atomic cloud with a central 
wavevector fe, and let us assume the polarization is one of 
the vectors V±. Then the maximum of the wave packet 
propagates with the group velocity Vg = Rc(Vfca;(fc)). 
As uj{k) depends on the angle 9 between B and k, Vg pos- 
sesses a component orthogonal to fc: in general, the wave 
packet does not propagate parallel to the wavevector. 
Since each polarization eigen-mode has its own direction 
of propagation, a birefringence effect takes place, as is 
well known in anisotropic media [35 . This magnetically- 
induced birefringence has already been observed [36], 
though under conditions differing from the ones described 
in the present article. 

We have checked that, provided |i5| is not much larger 
than ^B, the deviation of the wave packet from the direc- 
tion of k is negligible: the walk-off angle remains smaller 
than 0.05/(fc£o) [37]. Birefringence effects will thus be 
neglected in the following. 



E. Propagation in real space 

To study the propagation properties of the coherent 
mode in the atomic cloud, we need the Green's function 
of light in real space, Guj{r). It is the Fourier transform 
of ^;(fc,a;) Q: 



GUr) = 



Sk 

(2^ 



Afc 



S(/c,fc) 



(38) 



Neglecting birefringence effects, the angular part of the 
integral can be calculated with a stationary phase ap- 
proximation around the direction f: 

G^{r)^ -^ArUr)Ar (39) 

k sin(fcr) dk 



r/2 

2£res('5+ir/2) 



ArT(PP')A, 



■(40) 



In the presence of a magnetic field, Nj. is represented by a 
matrix with an antisymmetric part proportional to t], as 
can be seen from (31 1 and (16). This implies that Guj{r) 



is not a symmetric operator: the transposition operation 
is equivalent to flip the sign of ?], or equivalently to flip 
the sign of B: 



GUr,B)^'G^{r,-B) 



(43) 



This property is closely related to the reciprocity theorem 
|38| . As a simple illustration, assume a light beam is 
going successively through a linear polarizer e, the atomic 
cloud and a linear analyzer e'. The amplitude of the 
transmitted light is then proportional to Adir = e' • G^j ■ e. 
If we now consider the reverse situation where the light 
beam is traveling through the system along the opposite 
direction, its transmitted amplitude will now be A^ev — 
e • • e' = e' • *G(^ • e. We see that Adir = A^ev if 
and only if G^{r) is symmetric, which is not the case in 
the presence of the magnetic field. This is in essence the 
underlying principle behind optical isolators (or optical 
"diodes"), which are devices realizing A^ev — while 
Ad„ ^ 0. 



F. Magneto-optical effects 

The impact of B on the coherent propagation of light 
is embodied in the refractive index tensor ( 42 ) 



already mentioned in III C 



As 



it concerns a differential de- 
phasing and attenuation of the eigenmodes of propaga- 
tion. These magneto-optical effects have been extensively 
studied when light propagates parallel (Faraday effect) or 
perpendicular to (Cotton-Mouton effect) B. We present 
below both effects in terms of our formalism, giving sub- 
stance to the physical interpretation of the parameters C, 
•q and ^. 
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the traveled distance L. Hence, the parameter rj describes 
the Faraday effect. As one can also see, the parameter 
C plays the same role for the two eigenpolarizations : it 
acts as an isotropic refractive index. 

At low magnetic fields {iiB ^ F), 77 is proportional to 
B, and so is the rotation angle Q. The proportionality 
constant between this angle and the product BL is known 
as the Verdet constant Vb- For the F = 3 ^ Fe = 4 
transition of ^'^Rb, one finds 



Vb = 



3 

'4f4 



-8.10^ rad/(Tm) 



(45) 




B 

1 



(b) 




FIG. 1. (a) Resonant Faraday effect (5 = 0). A linearly- 
polarized light beam propagates parallel to B. The polariza- 
tion eigenmodes are the circular polarizations. In the course 
of propagation the polarization of the light beam keeps lin- 
ear but rotates around B (b) Resonant Cotton-Mouton effect 
(5 = 0). A linearly-polarized light beam propagates perpen- 
dicular to B. The polarization eigenmodes are linear, one 
being parallel to B and the other one being perpendicular to 
B. In the course of propagation the polarization of the light 
beam keeps linear but rotates around the propagation axis. 



1. Faraday effect 

The Faraday effect occurs when a light beam with a lin- 
ear polarization propagates along B in the atomic cloud. 
For sake of simplicity, we consider here that the frequency 
of light is exactly at resonance with the atomic transition 
{S — 0) so that 4> = 4>B, C, 1] and ^ are all real in the fol- 
lowing discussion. From ( 32 1 and ( 33 1 , the eigenmodes 



and their associated eigenvalues are: 



V 



A±^C±i\v\ 



(44) 



The eigenmodes thus identify with the left and right 
circular polarization vectors (see Figjl^). Their index 
mismatch ANr — ~\ri\/{k£o) is real, meaning that the 
two eigenmodes develop a phaseshift in the course of 
propagation. An elementary calculation shows that the 
polarization of the traveling beam remains linear but ro- 
tates around B by an angle 8 = riL/{2£o) proportional to 



with Iq ^ 100/im. This value is three orders of magnitude 
larger than those of classical materials, and comparable 
to the value measured in 1231. 



2. Cotton-Mouton or Voigt effect 

This effect describes the modification of the polariza- 
tion of a light beam propagating perpendicularly to B. 
As for the Faraday effect, we assume 5 = to simplify 
the discussion. The eigenmodes and their eigenvalues 
now read: 



V+ = z A+ = C + C V =y A_=C 



(46) 



The eigenmodes are now the linear polarization vectors 
(see FigjTjD). Their index mismatch AAV — iC/(2fc^o) is 
now purely imaginary, showing that V+ gets more atten- 
uated than V- by a factor e"^^/^^^"^ after traveling the 
distance L. Here again, the parameter C, plays the same 
role for the two propagation eigenmodes and acts as an 
isotropic refractive index. 

If the incident polarization is linear, it can be written 
as a real linear combination of and VL. During the 
propagation, the component along V+, i.e. along iJ, will 
decrease more than the one along i. e. perpendicular 
to B. As a whole, the polarization of the light beam 
remains linear, but rotates around fc by an angle which 
depends on 5, which is thus the parameter describing the 
Cotton-Mouton effect. 

In ref. [3S], the Cotton-Mouton effect is rather de- 
scribed as the transformation of an incident linear po- 
larization into an elliptical one in the course of propa- 
gation, this transformation being a consequence of the 
accumulated phase shift between the two propagating 
eigenmodes. This apparent contradiction can be lifted if 
one realizes that, contrary to our discussion, the Cotton- 
Mouton effect described in [3S] is in fact the one usually 
observed for light frequencies which are far-detuned from 
any resonance frequencies. Work out our theory at a very 
large detuning 5, we indeed recover the description given 
in [SSj. 
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3. General case 

When the direction of propagation k is neither along 
nor perpendicular to B, both effects mix. The Faraday 
effect will dominate when k is roughly parallel or anti- 
parallel to B, i. e. essentially when k is well inside the 
cone with apex angle for which (37) vanishes. On the 



contrary, the Cotton-Mouton effect dominates when k is 
essentially well outside this cone, i. e. roughly perpen- 
dicular to B. To give orders of magnitude, rj and ^ are 
comparable when </) ~ 1, which corresponds to i? ~ 2G 
for ^^Rb atoms. The apex angle is then 9q « 65°. The 
Cotton-Mouton effect then dominates for directions of 
propagation making an angle between 65° and 90° with 
B. In classical media, this would happen only in a nar- 
row angular width of order 10~^rad around the direc- 
tion orthogonal to B. This big difference in orders of 
magnitude is due to the strong resonant character of the 
atoms. All in all, the giant Faraday effect and the large 
zone of preponderance of the Cotton-Mouton effect make 
it necessary to take both effects into full account when 
studying the coherent propagation of light. 

When (5 0, Cj ^nd ^ are complex valued. Then 
the index mismatch is neither purely imaginary nor real. 
The two polarization eigenmodes still experience different 
phase shifts and different extinctions but the calculations 
and physical pictures lack the previous simplicity. 



IV. THE CBS EFFECT 

A. Independent scattering approximation 

At low optical density nX'^ <C 1, a semi-classical de- 
scription of propagation along scattering paths consist- 
ing of rays between consecutive scatterers is justified. 
For resonant scatterers, it implies fc^ 3> 1. As a con- 
sequence, scattering paths involving different scatterers 
are uncorrelated and the associated interference averages 
to zero. Recurrent scattering sequences (visiting a given 
scatterer more than once) can also be neglected defining 
the independent scattering approximation (ISA) (341. In 
this regime, the wave amplitude A is constructed as the 
coherent superposition A — A-p of the partial waves 
scattered along all quasi-classical scattering paths V join- 
ing the positions of the scatterers. Between two succes- 
sive scatterers the partial waves experiences the effective 
medium. In the ISA regime, the scattering amplitude A-p 
is thus computed using two building blocks, the scatter- 
ing by an individual atom and the coherent propagation. 



B. Multiple scattering and interference 



The average intensity of the wave / = (IX^p^'Pp) 
breaks into an incoherent contribution li = X^-pd^'Pp) 
and a coherent contribution Ic — 2Re(^p (A-p/ Ap)). 



The incoherent contribution itself breaks into the sum of 
a single scattering contribution Is and a diffuse one, Id, 
involving scattering paths containing more than 2 scat- 
terers. All these contributions depend on the polarization 
of the incoming light and on the detected polarization of 
the outgoing light. 

As is well known, the disorder average does not scram- 
ble two- wave interference effects between scattering loops 
traveled in opposite directions O EHl HO] ■ This is at the 
core of the CBS effect where interference between ampli- 
tudes associated to reverse scattering paths V and V (i. 
e. paths with the same sequence of scatterers but trav- 
eled in opposite order) contribute a constructive interfer- 
ence in a narrow angular cone around the backscattering 
direction [51|31I101- 



C. Backscattered intensity and the CBS contrast 

In this section, we compute the amplitude of multiple 
scattering paths and the backscattered intensity. We use 
the contrast of the interferences to determine the degree 
of coherence of light, which we express in term of a phase 
coherence length. 

Denoting by i? the angle between fe and the outgoing 
wave vector A;', the total average backscattered signal is 
I{d) = IsW + IdW + IcW where: 



I^{d) = 2 ^ Rela^a^ 



(47) 



V>2 



Here V > 2 means that scattering paths with at least 
two scatterers are included in the sum, Rp being the 
vector joining the endpoints of path V. The CBS signal 
Ic varies on a very small angular scale ~ ^/{ki) ^ 1, 
whereas the angular variations of Is and Id follows the 
Lambert's law and takes place over an angular range of 
order 1 radian. The incoherent contributions appear to 
be constant at the angular scale of the CBS cone and can 
be safely evaluated aX d = 0. 

As a two- wave interference, the CBS signal gives access 
to the degree of coherence of the outgoing wave and, in 
turn, at the coherence length of the scattering medium. 
The interference contrast is quantified by the CBS en- 
hancement factor a = l + Ic/{Is+Id) computed at = 0. 
which is the ratio of the total intensity at exact backscat- 
tering to the total intensity out of the backscattering 
cone. As single scattering events do not participate to the 
interference process, they decrease the contrast even if no 
dephasing mechanism is at work. When Is can be made 
to vanish or negligible, the coherence loss is directly asso- 
ciated to the ratio Id Id- For classical scatterers. Is — 
in the helicity-preserving polarization channel and reci- 
procity arguments show that Ic = Id and a = 2 in the 
absence of a magnetic field [131 UHl El] • When a magnetic 
field is present, Ic < Id and the coherence length of the 
medium becomes finite [20H22] . The situation for atoms 
with nonzero spin in the ground state proves more subtle 
and will be addressed in the next Section. 
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D. The intricacies of scattering under a magnetic 

field 

Consider a scattering path V containing s atoms, lo- 
cated at Ti, whose initial and final magnetic numbers 
are rrii and m[ {i = l,...,s). The incoming light angu- 
lar frequency is uj and its polarization vector e. Along 
path V, light propagates to the first atom, is scattered, 
propagates to the second atom and so on. After the s-th 
scattering event, light exits the medium and is detected 
in the polarization channel e'. One crucial point is that 
the internal Zeeman state of an atom with a degenerate 
ground state can change under scattering, which means, 
when a magnetic field is present, that the scattered pho- 
ton can have a different frequency that the incoming one. 
In other words, single scattering under a magnetic field 
is inelastic, the frequency change being of the order of 
liB. When /iS is larger than or comparable to S and/or 
r, this effect is not negligible: the effective medium and 
the magneto-optical effects felt by the photon depend on 
its frequency. 

If uj' denotes the outgoing frequency, the amplitude 
associated to path V reads: 

A-p ^e' ■ G^'{r' ~ rs)tm'^m,G^^_^^{rs - r^-i) 

• ■■G^^ ^{r2 - ri)t^'^^^G^{ri - r) • e (48) 

where 

i 

u)i^^+i^uj + g^iB'^{nia~m'^) (49) 

o=l 

is the frequency of light between the i-th and the {i + \)- 
th scatterers. The tensor tm'm has been defined in ([9|. 

For the reverse path V , the incoming angular frequency 
is still LJ and, by energy conservation, the outgoing an- 
gular frequency is still w'. The amplitude associated to 
V reads: 

Ap = e' ■ G^'{r - ri)t^,^„,^G^^ ^{ri - rs) . . . 

G^,,,_i(n-i - rs)tm'^m.,G^{rs - r') • e (50) 

with 

p 

=uj + giisB ^ {ma ~ m'^) (51) 

a—i+l 

The angular frequency of light traveling between atoms 
i and (i + 1) is Wi^i+i for path V and LOi+i^i for path V . In 
general, they differ and satisfy uji^i+i + w^+i^i = w -f w'. 
As the magnetic field introduces an explicit difference 
between A-p and Ap, the interference contrast will be 
reduced unless the change of frequency does not occur or 
is unlikely. This happens for example when the atom does 
not change its internal state under scattering (m^ — m[). 
The conditions for this situation will be examined in the 
next Section. 

As a consequence of this frequency change, the average 
over the internal degrees of freedom involves the whole 



scattering path. Indeed, because of the magneto-optical 
effects and the frequency change under scattering, the 
atomic internal and external degrees of freedom are intri- 
cated in a complicated way: after a scattering event takes 
place, the location of the next one depends on the value 
of the scattering mean free path, hence on the frequency 
of the emitted photon, hence on the change or not of in- 
ternal state. In such a situation, /j. Id and Ic can only be 
computed numerically. In practice one calculates the di- 
mensionless bistatic coefficients 72, = AirV'^uj'^Ix/i^i^'^A) 
[x — s,d,c), where A is the illuminated area. 

E. Monte-Carlo simulation 

Analytical results about the properties of the CBS cone 
can be obtained only in specific cases. In the absence of 
a magnetic field, the problem has been exactly solved for 
vector waves in a random medium of Rayleigh scatterers 
with a uniform density and a slab geometry [JT]. Fol- 
lowing the same lines, the solution has been extended to 
quasi-resonant atomic scatterers with degenerate ground 
states [12] • For other geometries, numerical calculations 
are necessary 

In the presence of a magnetic field, reference [22] con- 
tains a generalization of the analytical methods devel- 
oped in [IT] for Rayleigh scatterers but unfortunately 
fails to describe some aspects of the experimental results 
reported in [2(T . This has to be related to the approxi- 
mations done to compute the CBS cone, and in fine to 
the complexity of the exact calculation. For atomic scat- 
terers, the intrication between the external and internal 
degrees of freedom makes it almost impossible. 

The average intensity contributed by a 

scattering path V with s scatterers contains an average 
over the positions of the scatterers, i.e. a 3s-uple 
integral, and an average over the internal degrees of 
freedom, i.e. a sum over the final and initial Zeeman 
sub-levels of each scatterer in the ground state. To 
compute this multiple integral and these sums, we use 
a Monte-Carlo simulation able to extract at the same 
time Ap'' and A^' for a large number of paths. This 
numerical simulation allows us to take into account some 
experimental constraints, such as the shape and the 
density profile of the atomic cloud or the finite spectral 
width of the laser probe. Its principle is as follows: 

Step 1. A photon with frequency w, wave vector k — ku 
and polarization e enters the atomic cloud. It propagates 
on a distance r along u chosen according to the proba- 
bility distribution 




Here £~^{r) — n{r)a{(t)) is the inverse local scattering 
mean free path. It depends on the position of the photon 
if the atomic number density n is non-uniform. The 
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propagator is the Green's function in the real space (41 1. 
Step 2. The photon at position ru from the entrance 
point in the cloud is scattered by an atom. The atomic 
initial and final states \Fm) and \Fm') are chosen 
randomly with a uniform probability distribution. The 
scattering operator is t,n'm, eq.([9]). It transforms the 
incident polarization into the scattered polarization. 
Step 3. The photon angular frequency is changed by 
gfiB{m — m!). 

Step 4. A contribution to the single scattering bistatic 
coefficient 7^ is computed: the scattered photon is 
propagated along the backscattering direction § until it 
exits the atomic gas yielding the amplitude from which 
a contribution to 7^ is obtained. 

Step 5. The scattered photon is propagated towards a 
second scatterer. The direction of propagation is chosen 
according to an isotropic probability distribution to 
save computation time. The propagation distance is 
computed with the help of the distribution law eq.(52), 



smaller than 5% for magnetic fields up to {fiB/T ~ 10) 



for the results presented in section V C 



the change of frequency and polarization being taken 
into account. 

Step 6. The doubly scattered photon is propagated 

along the backscattering direction until it exits the 

atomic gas yielding the amplitude from which a contri- 
(2) 

bution to 7^ is obtained. 

Step 7. A photon, identical to the incident one, enters 
the atomic cloud, propagates along the previous double 
scattering path in reverse order, and exits the medium in 
the backscattering direction iJ. The scatterers experience 
exactly the same atomic transitions. This yields the 
amplitude associated to the reverse previous double 
scattering path from which, together with Step 6, a 
contribution to 7c is obtained. 

Step 8. The process is continued (triple scattering, etc) 
until the photon finally exits the atomic cloud. 
Step 9. Another incident photon is sent in the cloud 
and the whole process is repeated as many times as nec- 
essary to obtain a good signal-to-noise ratio. Typically, 
one needs to launch between 10^ and 10^ photons to 
obtain well converged values for 7^, 7^ — X^s'Td*^ ^'^^ 

7c = T,src '■ 



Up to the statistical errors, this method is quasi-exact 
and limited only by computer resources in the limit kl 3> 
1. 

When the magnetic field vanishes, the results of the 
Monte-Carlo simulations reported in 03] are recovered. 
At large magnetic fields, the modulus of the amplitudes 
associated to reverse paths are very sensitive to the scat- 
tering parameters of the paths. Any change in a direction 
of propagation modifies significantly the refractive index 
of the effective medium. As a consequence, the Monte- 
Carlo simulation needs to average over more and more 
fluctuating quantities when the magnetic field increases. 
The statistical error on the total diffuse intensity can be 
estimated by its standard deviation. It remains smaller 
than 1% for small magnetic fields {fiB/T < 1), and 



V. RESTORATION OF THE CBS CONTRAST 

In the following, we apply the results of the previous 
Sections to compute the CBS cone for quasi-resonant 
light propagating in a cold ^^Rb cloud. The frequency 
of light is chosen close to the frequency of the F = 3 — 
Ff,— A transition of the D2 line (wavelength A = 780nni, 
linewidth r/(27r) = 5.9MHz). The Lande factors of the 
ground and excited states are g = 1/3 and ge = 1/2. 
For this particular transition, a Zeeman shift = T 
corresponds to i? = 4.2G. 

At i? = 0, CBS experiments have reported very low en- 
hancement factors, e. g. a ~ 1.05 in the helicity preserv- 
ing channel |44[ I45j . This is in marked contrast with ex- 
periments with spherically-symmetric classical scatterers 
where reciprocity guarantees a takes its maximal value 2 
in the same polarization channel [IBj. A detailed analy- 
sis shows that the low a value observed with cold atoms 
comes from an imbalance between the amplitudes asso- 
ciated to reverse paths [T3] . This imbalance is noticeably 
absent for a F = — >■ Fe = 1 transition where a = 2 is 
recovered [47]. It is our goal in this section to show that 
the interference contrast can be fully restored with the 
help of an external magnetic field. 



A. Filtering out a closed transition 

The key idea to restore the CBS contrast in the 
helicity-preserving channel is simply to lift the degener- 
acy of the atomic ground state and to filter out a closed 
transition. This is done by splitting the Zeeman sub- 
levels with an external magnetic field (Zeeman effect) 
and by shining the atomic cloud with a light wave which 
is resonant with the |F = 3, to = 3) — |Fe = 4, me = 4) 
transition. To achieve this, one needs to impose 5 = 
(4^6 — 'ig)lJ'B = fiB. This transition is closed since an 
atom in the excited state |44) can only make a transition 
to the ground state |33). At sufficiently large i?, the other 
Zeeman sub-levels of the ground and excited states are 
sufficiently split away and are out of resonance, meaning 
that the |33) — )■ |44) transition is isolated. Thus, at large 
-B, the atomic cloud consists of (i) atoms which are in 
the sub-state |33) and can scatter light, and (ii) atoms 
which are not in the sub-state |33) and cannot scatter 
light because the frequency is too far-detuned from the 
other transitions. These |33)-scatterers behave like effec- 
tive two-level atoms which can only absorb and emit ct+ 
radiation, i. e. light with positive helicity along B. 

Under these circumstances, it makes no difference for 
light to travel a scattering path in one direction or the 
other. The multiple scattering amplitudes associated to 
any path V and to its reverse partner V are equal and 
the CBS contrast is restored. This restoration is expected 
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to be most spectacular in the helicity-preserving channel, 
because it is in this channel that the contrast is the lowest 
without any magnetic field. If the incident light beam 
is parallel to B and B is sufficiently large, it is easy 
to see that one only gets a non-vanishing CBS signal 
in the helicity non-preserving polarization channel. We 
will thus choose in the following the " Cotton-Mouton" 
configuration where B is perpendicular to the incident 
light beam and analyze the CBS signal in the helicity- 
preserving polarization channel. 



B. Small magnetic fields 

From the previous discussion, it seems that the con- 
trast restoration only occurs at sufficiently large B. In 
fact, it turns out that the contrast restoration even starts 
at small magnetic fields and gets larger as B is increased. 
To demonstrate this, we study analytically the single 
and double scattering signal in a uniform, semi-infinite 
medium in the limit jiB <^ T (meaning B 2.1 G for 
^^Rb). It is then possible to neglect the magneto-optical 
effects, and to propagate photons with the propagator 
(41) evaluated at i? = 0. This will be justified be- 
low in Section |VD| Expanding the scattering matrices 
Un'm at second order in (f)B — 2/ii3/r, and noticing that 
at same order, we then use expressions 



( 48 1 and ( 50 ) to compute the single and double scatter- 



ing amplitudes amplitudes. In the chosen geometry, the 
average over the external degrees of freedom when com- 

("2) (2) I I 

puting 7s, 7^ and 7c can be done analytically [28] . 
The calculations have been made here with the symbolic 
calculation software Maple^^ and yields: 



0.305 -fO 
0.217 0.14301 
1.166 + 0.04801 



(53) 
(54) 
(55) 



As one can see, both quantities increase with _B, meaning 
that the coherence length of the system is increased. 



C. Monte-Carlo simulations 

When the magnetic field is neither small nor large, 
there is no simple approximation that allows to compute 
analytically the bistatic coefficients, but they can at least 
be computed numerically with the help of the Monte- 
Carlo simulation described in Section IIVEI It is then 
possible to take into account a more realistic model of the 
atomic cloud than a semi-infinite uniform medium. We 
present results here for a spherically-symmetric atomic 
cloud with gaussian density and optical thickness 6 = 31 
(measured at B = and 5 = 0). We take a laser probe 
beam with spectral width equal to 0.3r. This allows a 
realistic comparison with our experiment done at a fixed 
total number of atoms. The computed values of a are 
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FIG. 2. Plot of the CBS enhancement factor a as measured 
in the parallel helicity channel h\\h for different values of B 
(circles) for light backscattered by a cold *^Rb atomic cloud 
in the Cotton-Mouton configuration fc _L B. For each value 
of B, the light is tuned on resonance with the |33) — >■ |44) 
transition [5 = /i-B). The spherically-symmetric atomic cloud 
is characterized by a gaussian density and an optical thick- 
ness b = 31 when 5 = and 5 = 0. One witnesses a dramatic 
increase of the CBS contrast compared to the situtation at 
5 = despite the fact that the time-reversal symmetry is 
broken. The solid line is the result of the Monte-Carlo simu- 
lation with no adjustable parameters. 



compared to the experimental ones for various values of 
B in Figj2j As one can see, a increases with B, starting 
from a = 1.05 at B = up to a ~ 1.35 at B = 40 G. 
The agreement between theory (solid line) and experi- 
ment (circles) is quite satisfactory. This shows that the 
Monte-Carlo simulation contains the essential ingredients 
that play a role in the restoration of the contrast. In the 
following, we will rest on the results of the Monte-Carlo 
simulation to elucidate the mechanisms at work by com- 
puting quantities which are not accessible to experiment, 
e. g. the bistatic coefficients for each scattering order. 

In Fig[3^, we plot 7c/7d as a function of B. This ratio 
is a measure of the degree of coherence of the outgoing 
light. This ratio grows when B increases, and tends to 1 
for large B (not shown in the figure). This confirms that 
the contrast of the interference, and in turn the coherence 
of the outgoing light, is actually restored by a magnetic 
field. 

In FigjSp we show how 7^ (triangles) , 7^ (crosses) and 
7c (circles) change with B. strongly decreases strongly 
with B because the atomic scattering cross section itself 
decreases. Meanwhile, 7c increases for magnetic fields up 
to 8 G. This shows the efficiency of the mechanism restor- 
ing interference. At larger fields however, the decrease of 
the scattering cross section takes over and 7c decreases 
again, although slower than 7^. At large B, 7c and 7^ 
tends to the same value (yielding a perfect coherence) 
but are outgrown by 7s, meaning a < 2. This behavior 
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FIG. 3. Results of the Monte-Carlo simulations under the 
same experimental conditions as Fig[2] (a) The coherence ra- 
tio 7c/7d increases with B. (b) Plots of the bistatc coefficients 
7s (triangles), 7d (crosses) and 7c (circles) in the backward 
direction as a function of B. The solid lines are drawn to 
guide the eyes. As one can see, 7c first increases than de- 
creases with B without varying too much. At the same time, 
7tj decreases strongly. At large B, jc = Jd and the coherence 
is restored. However the enhancement factor a < 2 because 
7s is not negligible. 



is not generic: at larger optical thickness, 7^ would have 
been smaller than 7^;. 

We also mention that the values of the bistatic coeffi- 
cients are independant of the value oi k£ ^ 1. However, 
the angular width of the backscattering cone depends on 
M. 



D. Impact of Faraday and Cotton- Mouton effects 

To study the impact of magneto-optical effects on the 
enhancement factor, the simplest way is to discard them 
in the Monte-Carlo simulation, and compare the obtained 
result with the experimental data in Fig[2] The magneto- 
optical effects distort the polarization of a propagating 
wave. Discarding them means here that a polarization 
is propagated without distorsion, only with attenuation. 
To ensure energy conservation, the attenuation length 
(i.e. the scattering mean free path) must be equal to £ — 
l/{n{r) a {(/))) where n{r) is the local density of atoms and 
cr(0) the total scattering cross section of an atom given 
by equation ( [l9| . All other parameters in the simulation 
are left unchanged. 

Fig|4] shows the plot of this expurgated enhancement 
factor a as a function of B (dotted line), together with 
its quantitative comparison to the real theoretical curve 
borrowed from Fig[2](sohd line). For B < AG {fiB < F), 
the impact of Faraday and Cotton-Mouton effects is neg- 
ligible. For larger fields, the true value a is slightly lower 
than a, the two curves being roughly parallel to each 
other. This means that magneto-optical effects do de- 
crease the phase coherence of the sample but this detri- 
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FIG. 4. Impact of magneto-optical effects occurring dur- 
ing propagation (Faraday and Cotton-Mouton effects) on the 
CBS enhancement factor a. The solid line is the theoretical 
curve obtained in Fig[2] The dotted line is the enhancement 
factor a calculated by discarding the magneto-optical effects, 
as one can see, the Faraday and Cotton-Mouton effects do 
decrease the contrast but their detrimental effect is counter- 
balanced and beaten by an efficient mechanism restoring the 
contrast. This mechanism is the modification of the scatter- 
ing properties of the atoms which, because of the Zeeman 
splitting, behave more and more as effective two-level atoms 
when B is increased and the light is tuned on resonance with 
a closed transition. 

mental effect is counter-balanced and beaten by a more 
powerful mechanism restoring coherence. As a matter of 
fact, at large B and after the first scattering event, a sin- 
gle eigenmode can propagate in the scattering medium, 
exemplifying why phase or extinction differences between 
the polarization eigenmodes cannot scramble the con- 
trast. Thus, the phenomenon explaining the restoration 
of the CBS contrast with B is really the modification of 
the scattering properties of the atoms which, because of 
the Zeeman splitting, behave more and more as effective 
two-level atoms when B is increased and the light is tuned 
on resonance with a closed transition. This is in marked 
contrast with classical scatterers where no such mecha- 
nism counter-balancing the detrimental magneto-optical 
effects does exist. 



E. Influence of higher and higher scattering orders 

When 5 = 0, the CBS effect observed with atoms 
having a degenerate ground state is dominated by dou- 
ble scattering paths, while higher-order scattering paths 
contribute significantly to the diffuse background [IS]. 
This is shown in the first line of Table |I] As B increases, 
the contrast is restored and higher and higher scattering 
orders contribute significantly both to 7c and 7^, see sec- 
ond line of Table [T] However, at the same time, the opti- 
cal thickness of the atomic cloud decreases. Higher-order 
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Scattering order 


2 3 4 5 


B=30G 


0.21 0.08 0.04 0.02 
0.81 0.75 0.69 0.38 



TABLE 1. CBS coherence factor ■yd 'yd for the first scattering 
orders when B = OG and B = 30G. The incoming hght and 
atomic cloud parameters are given in Fig[2] The atomic cloud 
is spherically symmetric with gaussian density. 




20 

B (Gauss) 



FIG. 5. Influence of higher and higher scattering orders 
on the CBS enhancement factor. The incoming light and 
atomic cloud parameters are given in Fig[2] We compare the 
theoretical enhancement factor borrowed from Fig|2| with the 
one obtained by only considering single and double scattering 
orders. At small B, high scattering orders contribute mostly 
to the diffuse signal. At large B, they contribute equally to 
the coherent and diffuse signals, making the CBS cone height 
increase. 



scattering paths become less and less probable and the 
two effects compete. The dotted line in figure [5] shows 
the enhancement factor calculated from the single and 
double scattering contributions alone. This approxima- 
tion overestimates the height of the backscattering cone 
at small magnetic fields, but underestimates it at large 
magnetic fields. This shows that high scattering orders 
do contribute to the CBS cone and cannot be discarded 
for a quantitative comparison. 



F. Role of optical pumping 



In section al D[ we mentioned that our theory does not 
take optical pumping into account (though it could be ex- 
tended to do so). In the present section, we give conclu- 
sive evidence that optical pumping is indeed negligible in 
our experiment by measuring the coherent transmission 
of the atomic cloud. The results are presented in Fig|6j 
for an incoming wavevector perpendicular to B and a cir- 
cular incoming polarization. If J = 0, the coherent trans- 
mission varies in time and its stationary value is shown in 
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FIG. 6. Coherent transmission of the atomic cloud as a func- 
tion of B for an initial optical thickness 6 = 31 measured at 
5 = and 5 — 0. Crosses: experiment for 5 = 0. Dashed 
line: theory for 5 = 0. Solid circles: experiment for 5 = /iB. 
Solid line: theory for 5 = ^B. When 5 = jj,B the incoming 
light is kept on resonance with the closed atomic transition 
|33) — >■ |44). When 5 = 0, the coherent transmission increases 
in time and converges to the plotted value. In this case, opti- 
cal pumping is at work, a situation not accounted for by our 
theory. For 5 = iJ,B, no time variation of the coherent trans- 
mission is observed. This is a strong indication that optical 
pumping is negligible in this case. 



Figj6]as a function of B (crosses). Our Monte Carlo sim- 
ulation (dotted line) is unable to reproduce these results 
for i? > 10 G, indicating that optical pumping is indeed 
present in our sample when 6 = 0. However, when the 
incident light beam is kept at resonance with the atomic 
transition |33) — >■ |44) (i.e. 6 — fiB), no time evolution of 
the coherent transmission is observed. This shows that 
the populations of the various Zeeman substates are al- 
most constant. The experimental data (circles) are well 
reproduced by the Monte-Carlo simulation (solid line). 
This is a strong indication that optical pumping is in- 
deed negligible when the incoming light is continuously 
kept at resonance with the closed atomic transition. 



G. Coherence length 

The notion of phase coherence length is a very im- 
portant concept in mesoscopic physics. It is the length 
scale at which, because of some dcphasing mechanisms, 
the interference effects as produced by the medium are 
effective. The larger is L^, the stronger is the impact of 
interference, and, in the case of the CBS effect, the larger 
is the CBS contrast. In the case of cold atoms, at B = 0, 
the degeneracy of the atomic ground state causes a loss of 
phase coherence between reversed scattering paths giving 
rise to a finite value of of the order of few mean free 
paths £ |13j. The increase of the CBS contrast when a 
magnetic field is applied is accordingly accompanied by 
a growth of L^. The Monte-Carlo simulation allows us 



15 



15 




I ■ ' ■ ■ ' ' ' ^ 

10 20 30 40 

B (Gauss) 



FIG. 7. Plot of Scjy — L^/l as a function of B, as extracted 
from our Monte-Carlo calculation. The incoming light and 
atomic cloud parameters are given in Fig|2] The points dis- 
persion reflects the numerical accuracy. As one can see, 
increases with B, roughly linearly (the solid line is drawn to 
guide the eye) . This is due to the lifting of the Zeeman degen- 
eracy which make atoms behave like effective two-level atoms 
when the light is tuned on resonance with a closed transition. 
This behavior is in sharp contrast with classical scatterers 
where s<^ decreases when B increases. 



to estimate the phase coherence length in the following 
way. In the presence of dephasing, the interference term 
7c associated to two reversed scattering paths of length 
L is related to the diffuse term 7^ by: 

7c~7de-^/^* =7de-^/^^ (56) 

where L/t = s, s being the scattering order, and L^/£ = 
S0. In Fig[7]we plot function of B as obtained 

numerically. It increases roughly linearly. 

One should note that our definition of the coherence 
length differs from the usual one where the distance trav- 
elled diffusively by the light inside the disordered sample 
is introduced [3j. With this convention, L cx y/s and 



H. Analogy with paramagnetic impurities in 
solid-state physics 

The surprising fact that a magnetic field can restore 
weak localization effects under well chosen circumstances 
although it breaks time-reversal invariance is already 
known in solid-state physics [48l |49]. In this context, 
one considers the propagation of electrons inside a metal 
at low (but finite) temperature, containing paramagnetic 
impurities. Thermal fluctuations make the spin of these 
impurities fluctuate in time. The scattering of an electron 
by such a fluctuating impurity randomizes the electron 
spin and the weak localization correction to electronic 



transport are reduced. This is similar to the loss of con- 
trast due to the degeneracy of the atomic ground state. 
When a large enough magnetic field is applied to the 
metal, the spins of the impurities are all aligned along B. 
The fluctuations of the spin component of the electrons 
are suppressed and the weak localization corrections to 
transport are restored. 

Finally, in both cases, the magnetic field freezes the in- 
ternal degrees of freedom and restores interference effect. 
The main difference with our case is that in solid-state 
physics, the magnetic field populates a unique spin state, 
whereas the Zeeman substates of the atoms are equally 
populated. It could be possible to realize an atomic cloud 
almost containing only atoms in the ground state |33), by 
using optical pumping. However, this would only restore 
the interference between reverse double scattering paths 
|19) . Nevertheless, it would be possible to first populate 
the 1 33) state, and then to apply an external magnetic 
field. This would increase the number of atoms partic- 
ipating to the scattering of light. Multiple scattering 
would then play a more important role and the enhance- 
ment factor would be larger than reported in the present 
article. 

VI. CONCLUSION 

To summarize, we have accurately described in this 
paper the propagation of light in cold atomic gases in 
the multiple scattering regime where ki ^ 1 and in the 
presence of an applied external magnetic field. In this 
regime, a semi-classical description is well suited and 
transport is described through a series of individual scat- 
tering events separated by coherent propagation in an 
effective medium. We have studied the impact of the 
magnetic field on the scattering of light by atoms with 
a degenerate ground state (differential and total cross 
sections) and the magneto-optical effects (Faraday and 
Cotton-Mouton effects) embodied in the refractive index 
tensor of the effective medium. Our results generalize 
previous works [531 HI] ■ We have then applied our the- 
ory to the study of the coherent backscattering effect 
and we have shown that the magnetic field can lead to a 
full restoration of the two- wave interference contrast pro- 
vided the incoming light is continuously set on resonance 
with a closed atomic transition as B is increased. The 
reason for the restoration of contrast is that the atoms 
behave as effective two-level atoms for which scattering 
amplitudes associated to reverse scattering paths have 
the same strength. 
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